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Abstract 

Unidirectionally coupled dynamical system is studied by focusing 
on the input (or boundary) dependence. Due to convective instabil- 
ity, noise at an up-flow is spatially amplified to form an oscillation. 
The response, given by the down-flow dynamics, shows both analogue 
and digital changes, where the former is represented by oscillation fre- 
quency and the latter by different type of dynamics. The underlying 
universal mechanism for these changes is clarified by the spatial change 
of the co-moving Lyapunov exponent, with which the condition for the 
input dependence is formulated. The mechanism has a remarkable 
dependence on the noise strength, and works only within its medium 
range. Relevance of our mechanism to intra-cellular signal dynamics 
is discussed, by making our dynamics correspond to the auto-catalytic 
biochemical reaction for the chemical concentration, and the input to 
the external signal, and the noise to the concentration fluctuation of 
chemicals. 

1 Introduction 

In a biological system, signaling phenomena are important, that transform 
external inputs to outputs. In general the signaling process is not simple, 
and is not given by a fixed one-to-one correspondence. A typical example 
is an intra-cellular signaling, where several chemicals are involved. These 
chemicals include several catalytic reactions and positive feedback process 
to amplify the input. As a result nonlinear dynamics of coupled elements is 
concerned. 

In a cell, an external stimulus coming from the outside of a cell is selected, 
amplified, and transported gnal. Recent development in molecular 

biology provides us more and more detailed information on the signaling 
process[|I|. In the experiment of molecular biology, mechanisms on chemical 
signal transmission have become clearer. Network of the pathways from a 
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given stimulus through signal chemicals to the final response is also being 
clarified. 

The network of signaling pathway is generally very much complicated. A 
large number of biochemical reactions are involved, from the cell membrane 
to ligand, and to the final response reaction. Often, the pathway may not 
be uni-directional, but several pathways interfere with each other. With this 
complication one could even expect that there would be no correlation be- 
tween input and final response. If such were the case, however, no biological 
function would be expected. Then, why is the network of signal pathway 
so long and complicated? How can the signal pathway respond suitably to 
inputs? How does the signal process work within possible thermodynamic 
fiuctuations? 

To answer such general questions, experimental molecular biology is not 
sufficient, but some theoretical study on general features of a biosystem is 
required. To address the above questions, we start from a study of simple 
model pathway formed uni-directionally, as an ideal simplification of the 
signal pathway, where a model with a sequence of biochemical reactions from 
an input is adopted, that transforms input to the final response. Here we 
do not take a biologically realistic model, but focus on a dynamical system 
aspect abstracted by a simple model. With this we will propose a mechanism 
how a signal is amplified and transmitted within fiuctuations. 

To understand the mechanism of signaling, we note that the signal system 
should satisfy at least the following three properties; 

property 1 : Amplification and transmission of input ; external signal is 
amplified and transmitted into an intra-cellular dynamics. 

property 2 : (Linear) stability : the system returns to the original state 
when the input is off. 

property 3 : Input dependence : response following the internal dynamics 
depends on the nature and strength of the input signal. 

In the present paper, we introduce an abstract model with a simple bio- 
chemical signal pathway, and show that its dynamics can satisfy the above 
three properties. A chain of coupled biochemical reactions represents the 
model. Indeed, in signal pathways, there are several chains of reactions, 
where by means of enzyme-catalyzed pathways, many biochemical reactions 
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are driven in one direction, through couphng to energetically favorable hy- 
drolysis of ATP to ADP and inorganic phosphate [0. 

In a dynamical system with one-way coupling, an important notion is 
convective instability, where small perturbation is amplified to down-flow, 
although the system is linearly stable and provides a stationary state without 
such perturbation. We will show that the properties 1 and 2 are satisfied 
with this convective instability, as will be explained in the next section. It 
will be also shown that some of the one-way coupled dynamical system with 
convective instability can have input dependence. 

Here we consider a network of catalytic reaction elements, each of which 
is represented by differential equations showing excitatory dynamics. Each 
element is assumed to be connected uni-directionally. In other words, each 
chemical at the {i — l)-th site has an influence at the next z-th level. To be 
specific, we choose the model, 

Xi = F{X\X'-^) (1) 

with X* = {x\ y^). 

In the present model, the input is represented just as a boundary condition 
X^ = const., while the response is given by the dynamics at the down-flow 
i ^ 1. We will find input dependence, or in other words boundary-condition 
sensitivity, and clarify its mechanism in the term of convective instability. 

Since a convectively unstable (CU) system is often sensitive to noise as 
described later, it is important to discuss how our signal transmission works 
in the presence of noise, due to chemical fluctuations. Indeed, it is shown 
that our mechanism works in the presence of noise, or rather it works best 
when the noise strength is of a medium size. 

Although our motivation is originated in signaling phenomena, our formu- 
lation of input-dependent djTiamics is generally applied to any uni-directional 
chain of dynamical systems. In this respect, the mechanism of analogue and 
digital input-dependence, as well as our formula for their condition, will have 
general applicability. It may include neural networks 0, optical networks 0, 
open fluid flow, and other chain network of dynamical systems. 

The present paper is organized as follows. In §2, noise amplification by 
convective instability will be reviewed. In §3, a simple model using coupled 
chemical reactions is introduced for biological signaling pathway. In §4, input 
dependence in our model is presented, while its mechanism will be explained 
in §5 in connection with convective instability and its spatial change due 
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to input signal. In §6, the relevance of the noise to the input dependence is 
described. It is shown that the input dependence is possible only within some 
range of noise amplitude. §7 is devoted to the discussion and the conclusion. 

2 Convective instability in a one-way coupled 
system 

2.1 convective instability 

One-way coupled system is introduced @, |^ as an abstract model for an 
open flow system, for example, in a fluid system. By spatial amplification of 
disturbance at upper flow, the dynamics increases its complexity, as it goes to 
down-flow]^, ^, ^. This change of dynamics is often triggered by convective 
instability. 

Convective instability expresses how some perturbation is amplified along 
a flow. If a system is 'convectively unstable' (CU), perturbation is spatially 
amplified and transmitted as in Fig.[I|. On the other hand, if a system is 
convectively stable (CS), perturbation at an upper flow is damped as it goes 
down-flow. Note that even if the perturbation is damped at each site (i.e., 
the system is linearly stable (LS) ), the system can be CU as shown in Fig.|l|. 



Convectively Stable ( CS ) Convectively Unstable ( CU ) 




time Hiile 



Figure 1: Schematic representation of the evolution of perturbation 

Convective instability is quantitatively characterized by co-moving Lya- 
punov exponent A^,, Lyapunov exponent observed from the inertial system 
moving with the velocity f For a given state, if maxv Xy is positive, 

the state is convectively unstable. The convective stability is supported by 
maxy At, < 0, distinguishable from linear stability, given by the condition 
Ao < 0. The co-moving Lyapunov exponent is usually applied to an attrac- 
tor, where chaos with convective instability is characterized by the positivity 
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of maXy \y. It is however applied to any state, to characterize its stabihty. 

Since A^, characterizes the amphfication of perturbation for the velocity v, 
the amplification per one lattice site is given by — . Hence the amplification 
rate per one site is given by the spatial Lyapunov exponent A'^[p!0[]; 

= max — (2) 

V V 

2.2 noise sustained structure in convectively unstable 
system 

For a system with convective instability, noise plays an important role. When 
a system is linearly stable but convectively unstable, applied noise is spatially 
amplified and transmitted from up-flow to down-flow, until spatiotemporal 
structure is generated at the down-flow|Q, 0. This structure is different 
from that for a system without noise. 



Mechanism of the structure formation is summarized as follows [|], [IT 
Assume that noise is added to LS but CU fixed point. Around the fixed 
point, the noise is spatially amplified and transmitted from up-flow to down- 
flow. The more it goes to down-flow, the larger oscillation can appear, until 
some stationary dynamics (such as periodic oscillation) is generated for i ^ 
1. As long as the noise is added at the most up-flow element {i = 0), 
the down-flow dynamics remains same. This noise-induced structure in a 
convectively unstable system is a general feature in a one-way coupled system 
and important for our model. Of course, if the system is CS around the fixed 
points at all elements, noise is spatially damped, and no oscillation exists at 
the down-flow. 

3 Model 

As a simple signaling model, we choose a one-way coupled differential equa- 
tions(OCDE)Q, where differential equation of each element can be regarded 
to express the dynamics of biochemical reaction. 



^OCDE is also studied by Aranson et.al||]. 
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3.1 dynamics of a single element 



As dynamics of a single element, we choose an auto-catalytic biochemical 
reaction system, consisting of activator x and inhibitor y (see Fig.^). In other 
words, there is a set of two chemical variables for chemical concentrations, 
represented by differential equations with two degrees of freedom. 




activate 



X : activator 
y : inhibitor 



K 

reservoir 



Figure 2: Schematic representation of biochemical reaction of a single ele- 
ment 



As a specific example, we choose the following model, 

x = f{x,y) 
y = 9ix,y) 

with 

f{x,y)=x{{l-x)ax-by} + K 
9{x,y) = y{cx -dy) + K 

where the parameters a, 6, c, d, and K can be related with the rate constants 
of biochemical reaction. In the present paper, we set K <^ 1, a,b,c,d ^ 1 
so that the system is satisfied with the three properties in §1. In the present 
case, the suppression of the catalytic process is provided by the term (1 — x) 
(where x is between and 1), but other forms like Michael- Menten's can be 
adopted. Here details of the model are not important, and only the excitatory 
nature in the dynamics is necessary. 

This single-element dynamics is chosen so that (x, y) converges to a lin- 
early stable fixed point (x*,?/*), where f{x^:,y^) = g{x^,y^:) = 0. There is 
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neither a stable limit cycle nor another fixed point besides the above fixed 
point. 

3.2 one-way coupling 

Assuming that the identical set of dynamics is set at each i-th node, the 
coupling from each node to the next node is introduced as the activation 
process of the i-th node by the activator of the (i — l)-th node. In other 
words, we assume that the activator chemical at the i-th node is catalyzed 
by the {i — l)-th activator. This leads to the following one-way coupled 
differential equations. 

xi =f{x\y\x'-^) 

=g{x\y') ^ ' 



f{x^,y\x^ ^) = x^ {{1 — x^){ax^ + ex^ ^) — by^} + K + rjl , . 
g{x\y') ^y'{cx'-dy') + K + rji 

where i denotes a spatial position or a level of signal pathway. Coupling 
constant e corresponds to the rate constant of biochemical reaction between 
(i — l)-th activator and i-th activator, and 77 represents white noise expressing 
concentration fiuctuation of chemicals satisfying 

{vi{tH{t-r))t = S,,,S,,kS{T)a' (7) 

with o", as the strength of the fluctuation, which is independent of i, x or 
y. All the parameters arc assumed to be independent of the node i, for 
simplicity. Again, the details of the coupling form are not important, as long 
as some nonlinear one-way coupling is included. 

This pathway transforms external signal into final response. Here, the 
input signal is given by the concentration x^ (external signaling chemical) 
that appears as a parameter for the dynamics of x^. This chemical concen- 
tration x^ is set as a constant. The response is given by the concentration 
of chemicals at the down-fiow (i 3> 1), whose dependence on x^ is studied. 
In other words, x'^ works as the boundary condition of the one-way coupling 
system. In the present paper, we do not consider its temporal change (like 
input oscillation), although its introduction will be an interesting problem in 
future. We study how the excitatory pulse appears at the down-flow. 
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With regards to the connection to cell signaling problem, the node i can 
be regarded either to one-dimensional spatial position or to a level of kinase 
reaction. In the latter case the use of identical reaction eq.(^ for each i is 
of course, too drastic simplification. Still, it is useful to discuss a possible 
signaling mechanism in this simple case. Indeed the mechanism and concepts 
in the present paper can be straightforwardly applied even if the reaction 
equation differs by each element i. 

3.3 numerical results : pulse generation 

Since we have chosen our model to satisfy the linear stability of fixed points, 
all elements converge to a fixed point pattern, when the noise is not added 
as shown in Fig.|^. The pattern shows dependence on the boundary at the 
up-fiow, but it rapidly converges to a fixed value. The down-flow pattern 
is spatially constant whose value does not depend on the input. There is 
no stable attractor with oscillatory dynamics besides the above fixed point 
pattern. 




Figure 3: Snapshot of space-amplitude plot for our system with noise (+ : 
cr = 10^^) and without noise (x : cr = 0). With noise, the dynamics of each 
element at the down-flow is periodic, a = 0.4, 6 = 5.12, c = 2.0, d = 3.55, 
K = 0.0004, e = 2.8. Throughout the figures in the paper these parameter 
values are fixed. 



Here we are interested in the case that the fixed point pattern is convec- 



tively unstable at least some point. Then the mechanism mentioned at § |2^ 
works when some noise is added. We have found that a pulse-like oscillation 
is generated at the down-flow, as shown in Fig.Q (see also Fig.^). Often, 
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(i.e., for most parameter values and input values), the generated oscillation 
is periodic. At the down-flow, the oscillation becomes periodic both in space 
and time. In other words, the pulse is transmitted to down-flow without 
changing its pattern (i.e., amplitude or frequency). Note that the limit cycle 
has to be CS to be transmitted to down-flow without being changed by noise 
added at each element. Still, this limit cycle disappears if the noise is turned 
off. 

4 Input dependence 

To see the input dependence of the pulse-like solution, we have changed the 
concentration x° and studied how the dynamics at the down-flow changes. 
Here the boundary (input) concentration is taken to be flxed in time, and 
temporal information of the input is discarded. 

Two types of the input-dependence are discovered, as to the down-flow 
dynamics, given by the behavior of {x^{t),y^{t)) at i ^ 1. The flrst one 
is a 'digital' (i.e., threshold-type) change, represented as a different phase of 
dynamics, such as flxed points, oscillation, and so forth, while the second one 
is 'analogue' change of the frequency or amplitude of the oscillation. Note 
that these input dependence of the down-flow dynamics are kept even for 
i — *■ cxo where the dynamics spatially converges. 

4.1 digital change 

With the application of noise to all elements, the response shows the fol- 
lowing three phases successively, with the increase of the input value (the 
boundary condition) The difference of the phases is shown in Fig.^ with 
the spatiotemporal pattern and the corresponding power spectrum for the 
time series of at the down-flow. 

periodic oscillation phase (p-phase) (limit cycle) 

At the down-flow, periodic oscillation with a large amplitude is gen- 
erated, as shown in Fig.^(a). The motion at the down-flow is quite 
regular, in the presence of noise, as is shown in the power spectrum of 
X* at 2 = 80 (see Fig.^(b)). The orbit of {x\y^) is plotted in the phase 
space in Fig.|^. Note that the orbit does not pass through the fixed 
point of the noiseless case. 



9 



As mentioned at if the noise is turned off, the oscillation damps 
and each element dynamics converges to a fixed point. On the other 
hand, the oscillation is CS and transmitted without influenced by noise. 



stochastic oscillation phase (s-phase) At the down-flow, stochastic os- 
cillation with a large amplitude is generated. The interval between the 
peaks is stochastic and longer than the period of oscillation in the p- 
phase (see Fig.H(c) for the spatiotemporal pattern). Amplification and 
transmission of noise are regular at the p-phase while they are possible 
only stochastically at the present phase. Peaks in the power spectrum 
of the down- flow dynamics are no longer observed, as shown in Fig.^(d), 
and are replaced by the broad band spectra. The orbit stays close to 
the same fixed point of the noiseless case for long time, emits a pulse 
train stochastically once, and returns to the former fixed point. (See 
Fig.^ where the orbit of (x*,?/*). ) 

Although the oscillation is aperiodic, each pulse itself remains to be CS 
and is transmitted without influenced by noise. 

fixed point phase (f-phase) No pulse with large amplitude is generated. 
The orbit stays around the fixed points with some fluctuation (~ a). 
No spatial amplification is observed. At the transition from the s-phase 
to the present phase, the frequency of the pulse generation goes to zero. 

We have measured the distribution of the interval between the peaks of 
X* at the down-flow, as the boundary x° is changed. At the p-phase, the 
distribution has a sharp peak, while no typical peaks are observed at the 
distributions for the s-phase. 

To see if the dynamics converges at a down-flow, we have measured the 
variance of x, that is < (x*)^ > — < x* >^, with < .. > as the temporal 
average. Fig.|^ gives the spatial change of the variance for three different 
boundary conditions x°. The upper plot corresponds to the p-phase, the 
middle to the s-phase, and the bottom corresponds to the f-phase. As in the 
figure, the dynamics spatially converges at i ~ 30, while the input depen- 
dence is preserved for i oo. (Indeed we have confirmed that the variance 
< (x*)^ > — < X* >^ remains different for each boundary condition.) All of 
our results show that the input dependence is not a transient in space but is 
kept for 2 — ^ oo. 
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Figure 4: Spatiotcmporal plot of x^(t) (a), (c) and the corresponding power 
spectrum (b), (d) with i = 80. The input value are 0.070 (a, b) and 0.054 
(c, d). (a) and (b) show periodic oscillation at down- flow, (c) and (d) give 
stochastic oscillation, cr = 1.0 x 10~^. Parameters a, 6, c, d, K, e are same as 
Fig.3. 




Figure 5: Phase space at i = 80 to different input value Parameters 
a, 6, c, d, e and intensity of noise a are same as Fig.4. The dotted line 
corresponds to Fig.4(c) (s- phase), while the sohd line corresponds to Fig.4(a) 
(p-phase). 
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Figure 6: Dependence of the variance < (x*)^ > — < >^ on the site i. 
The top one (o) corresponds to the input of Fig.4(a) (p-phase), the middle 
one (x) to Fig.4(c) (s-phase), while the bottom one (□) corresponds to the 
fixed point phase (f-phase). The dynamics spatially converges at i ~ 30 and 
input dependence is kept for i ^ oo. Parameters a,b,c, d, K, e and intensity 
of noise a are same as Fig. 4. 



4.2 analogue change 

Analogue dependence on input is given by a continuous change of the fre- 
quency and amphtude of the oscillation (at i ^ 1) with the input, i.e, the 
boundary condition It is observed both in the periodic and stochastic 
oscillation phases. Fig.|^ shows input dependence of the average frequency of 
the down-flow dynamics. In both the phases, the average frequency of the 
down-flow dynamics changes continuously with the boundary value x^. The 
down-flow dynamics shows not only digital but also analogue dependence 
on the boundary. We have also measured power spectrum of the down-flow 
dynamics. The position of the peak frequency (accordingly the response 
dynamics) continuously changes with the input x^. 



5 Mechanism 

How is the input dependence formed? Recall the following two properties 
described in § |2.2| : First, the formation of a pulse train structure depends on 
if the fixed point solution is CS or CU. Second, the nature of the generated 
oscillation is influenced by the degree of convective instability of the fixed 
point at the region where the pulse is being formed. 
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0.02 0.04 0.06 0.08 0.1 

Input (boundary condition x'^0) 



Figure 7: Input dependence of the average frequency of pulses at down-flow 
is plotted. From the power spectrum, it is shown that the p-phase lies at 
> 0.054, and s- phase at 0.054 > > 0.051. For both the phases, the 
average frequency changes with the input value x'^. Parameters a, b, c, d, K, e 
and intensity of noise a are same as Fig.4. 
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By referring to these properties, the mechanism of input dependent dy- 
namics is sketched as follows: Depending on the input, the fixed point values 
at the upper flow change, which leads to the difference in convective insta- 
bility. Accordingly the degree of the noise amplification changes, and the 
generated dynamics at the down-flow is different. In this section we study 
the above mechanism of the input dependence in the term of the convective 
instability of the fixed point. 

Without noise, all elements converge to fixed points {xl,yl) that depend 
on the site i. Now we introduce the site-dependent co-moving Lyapunov 
exponent \v{i) and spatial instability exponent X^{i) of {xl,yl). We will see 
that the input dependence can be explained by spatial change of for 
{xi,yi). 



5.1 spatial convergence of fixed point 

In our model, {xl,yl) is determined as follows. 

xi = f{xi,yl,xi-^) =0 

yl = 9{xi,yi} =0 

Eq.(||) gives a recursion equation of from yl~^) to (x* , yl). Given 
at the left boundary, this equation can be solved successively for i = 1, 2, .... 
by the iteration. The trace of spatial convergence, {xl,yl) — > {xl,yl) — *• 
{x'^,y'^) is decided as in Fig.^. In our model, {x'^,y'^) is decided uniquely^ 
to Vx°. We use the notation. 



(X 



*j y*) 



{xT,yT) 



= A^(oo) 



(9) 



Now let us study the convective stability of this fixed point. Since the 
fixed point is site dependent, we extend the co-moving Lyapunov exponent 
to a spatially local one, by measuring the amplification of perturbation to 
the next site from each site for a given velocity. Technically this amplifica- 
tion rate is measured as the co- moving Lyapunov exponent for a snapshot 
homogeneous pattern {x\y'') = (xl,yl) for / > i (see Appendix). 

Fig.|^(a) gives the co- moving Lyapunov exponent Xv{i) of the fixed point 
at each site, while of the fixed point at each site is given in Fig.^(b). 

^It is an interesting future problem to study a case with a multiple solution of {x^,y^). 
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Figure 8: The fixed point is plotted for different input values x^. Pa- 
rameters a, 6, c, d, K, e are same as Fig.4. Different marks corresponds to 
different input values, x° = 0.01(+), 0.02(x), 0.03(A), 0.04(n), 0.05(v), 
0.06(o), and 0.07(«). 




Figure 9: Co- moving Lyapunov exponents A^(i) (a) and the spatial instability 
exponents (b) are plotted as a function of for i = 0(0), 5(-|-), 10(n), 
15(x), 20(A), 25(*), 30(0), 35(+), 40(n). Parameters a, b, c, d, K, e are same 
as Fig.4. The boundary condition is given by ~ 0.060. is maximum 
at ~ 0.1 for every site. 
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Here, site dependence is clearly seen. At i < 30 the fixed point is CU, since 
Xv{i) is positive at f ~ 0.1. On the other hand, at i > 30, the fixed point 
is CS, Xv{i) is negative for all v. In our model, the fixed points change from 
convectively unstable to stable ones, as the site goes down-fiow. 

Recall the growth rate of the perturbation to the next sites is given by 
the spatial instability, X^{i) = maXy In our model, such velocity v that 
maximizes is independent of i, as is shown in Fig.|^(b). Change of X^{i) 
with the input is shown in Fig.[iy. In our model, the change is monotonic 
with respect to site i f\, either by decreasing from Af (1) to Af or increasing. 

Now we discuss the relationship between Af (i) and input dependence. In 
Table |I], we have classified our dynamics into three types according to the 
convective stability at the up-fiow and down-fiow. 

Table 1: Composition of CS or CU states at up-fiow and down-fiow 





fixed point 
at up-fiow 


fixed point 
at down-fiow 


limit cycle 
at down-fiow 


Type 1 


CS 


CS 


none 


Type 2 


CU 


CS 


CS 


Type 3 


CU or CS 


CU 


CS 



First, for the Typel, no oscillation is generated, since the fixed point is 
CS. The fixed point approaches a constant value for i ^ oo, and the state 
value at the down-fiow is input independent. For the Type3 the fixed point 
at the down-fiow is CU, and spatial amplification from a fixed point always 
occurs, independently of the input value Hence, no digital change occurs. 
Indeed, in numerical simulations, the digital change from a fixed point to 
oscillation is observed only at the Type2 case only. In this case, there are 
two CS states at i — oo. One is the fixed point {x^, y^) which is LS (Linearly 
Stable), and the other is a limit cycle. According to the nature of the up-fiow 
dynamics, either of the two states is selected, that leads to a digital change 
in the dynamics. 



■^It is interesting to study the case with non-monotonic convergence of fixed points 
(such as a spatial periodic or chaotic case), in future. 
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Figure 10: The spatial instability exponent Af(i) of the fixed point. Differ- 
ent marks represent different input value, = 0.01(+), 0.02(x), 0.03(A), 
0.04(n), 0.05(v), 0.06(o), and 0.07(«), corresponding to Fig.8. 
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5.2 mechanism of digital change 

In this section mechanism of digital change is quantitatively analyzed ac- 
cording to the CO- moving Lyapunov exponent. 

Since the fixed point becomes convectively stable at the down-flow, the 
oscillatory dynamics should be formed before the fixed point is stabilized. 
Hence the condition for the formation of oscillatory dynamics is given by the 
competition between two time scales, the scale iu where the lattice remains 
to be CU, and other the scale ig, required to the generation of a pulse, i.e. 
the "growth scale" for perturbation. The scale iu is given by the site where 
the fixed point (x* , yl) becomes CS, and it is given by the condition > 
at i < iu and X^{i) < at i > iu- 

The scale ig is estimated as the scale where the noise is amplified to the 
scale of 0(1), to generate oscillatory dynamics. As long as the noise is small, 
the amplification rate of noise is determined by the spatial instability Xf{i) 
of fixed point Hence ig is roughly estimated as follows. 



'9 

a ^ 

1=0 



ne^*» = ae&^*«~l (10) 



ixm-log- (11) 

Since the perturbation is no longer amplified at a CS state, it is necessary 
that noise is amplified while the dynamics at the site remains to be CU. The 
noise has to be amplified to 0(1) before i ~ iu, the condition ig ~ iu is 
imposed to have oscillatory dynamics. As the input (boundary) changes, 
fixed points (x*,?/*) change accordingly, and both i„ and ig change. Then 



the relation between ig and iu can change qualitatively. Fig.[TT| shows input 
dependence of ig and iu for our model for parameters belonging to Type2. 

Note that the relationship ig < iu changes with the input xq. For > 0.065, 
ig is smaller than iu, and indeed we have observed oscillatory dynamics with 
a large amplitude. Around ~ 0.065, ig ~ iu, where stochastic oscillation 
is found with intermittent pulse generation from the fixed point state. No 
pulse is generated for x^ < 0.065. 

We can summarize the change of dynamics according to the relationship 
between ig and iu, as follows: 
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• ig < iu '■ This gives periodic phase, where stationary oscillation is 
formed at down-flow. Since the generated oscillation itself is CS, the 
pulse is transmitted to down-flow without affected by noise. As long as 
noise is added stationarily (at least at the up-flow), spatially stationary 
state is formed at the down-flow. 

• ig ^ iu '■ This corresponds to the stochastic phase, where the for- 
mation of pulse is strongly sensitive to noise at the moment. Formed 
pulse itself is CS and can be transmitted to down-flow stationarily, but 
its formation is intermittent. Indeed, we have measured 4^ at the pa- 
rameters of the stochastic phase, and the value stays around 1.0 ~ 1.2. 

• ig = oo {ig > iu) '■ This give a fixed point phase, since the fixed 
point becomes CS before a pulse is generated. No pulse exists, and the 
down-flow dynamics remains around the fixed point with some noise. 
Indeed, the condition(^) is not satisfied for any ig. 

This summarizes how the digital change appears for our dynamics. It 
should be noted that the above mechanism generally holds in a one-way 
coupled system with the fixed points convectively unstable at the up-flow 
and stable at the down-flow. Now it is also clearer why the digital change 
is seen only for the Type2. For the Typel, iu is oo, and only the fixed 
point phase appears, while for the Type3, iu is zero, and the pulse is formed 
irrespective of the boundary condition (input). The mechanism is shown 
schematically in Fig.|12[ 

5.3 mechanism of analogue change 

Here we discuss the analogue change against input, that 'transforms' the 
input concentration to the response frequency. The rate of formation of 
pulse is expected to depend on the fixed points and the convective instability 
before the pulse train is formed. On the other hand, as the site goes down- 
flow, the fixed point approaches (a;*,?/^,) independent of the input In 
other words, the spatial instability is independent of the input, after the 
fixed point converges closely to (a;*,?/*). Hence, only if the scale ig, required 
to form the pulse, is smaller than the scale of the above convergence, the 
dynamics generated at the down- flow is expected to depend on the input. 
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0.03 0.035 0.04 0.045 0.05 0.055 0.06 0.065 0.07 0.075 0.08 
Input (boundary condition x°) 



Figure 11: Dependence of iu (+) and ig (o) on input value For the ref- 
erence the change of variance of x^^ (1^) is plotted. Parameters a, b, c, d, K, e 
are same as Fig.4. Intensity of noise a — 0.001. 



Digital Cliange 



ig< /„:p-phase 
ig~ /„:s-phase 
ig> /„:f-phase 




i„ change 



Figure 12: Schematic diagram of the change of Af , corresponding to each of 
the three phases. 
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As a simple measure of the relaxation length v of fixed points, we will 
introduce a 'half-decay' scale. Since the spatial convergence of Af (i) is mono- 
tonic as described in §5.1| and Af(z) characterizes the amplification of noise, 
we use the convergence of Af (i) to define ir- As the half-decay scale of Af(«), 
ir is defined such site that satisfies Af (v) = ^^i^^-^. Roughly speaking, 
Z)} Af («) changes sensitively on i ioi i < ir, while for i > ir, it weakly depends 
on i. This is expected, for example, by assuming that Af(i) exponentially 
relaxes to Af . 

Now let us focus on the relationship between ig and ir- If the convergence 
scale ir is much smaller than ig, the rate of amplification of noise given by 
J2) Af (j) depends little on the input. On the other hand, if v is larger than 
ig, the amplification rate largely depends on the input. Since this change 
of amplification rate is continuous with the input value, analogue change of 
response is expected. Hence the following two cases are classified. 

Case A : ig < ir 

At the region where the pulse is formed (i < ig), the spatial convergence 
of fixed points is not yet completed. The noise amplification rate there 
is dependent on the input. Thus the nature of generated wave ai i ^ ig 
depends on the input. Since the formed wave pattern is convectively 
stable, this dependence on the input is preserved to down-flow. 

Case B : > V 

Since the fixed point has almost converged to {x^,y^) at the site where 
the wave pattern is formed, the amplification of noise and the nature 
of formed wave are insensitive to the input value. Hence the nature of 
generated wave depends little on the input. 

To write down the condition for the analogue dependence for the Type2, 
we need to add the condition for the formation of wave (discussed in § p.2| ), 
while for the TypeS the wave is always formed. Hence the conditions of the 
analogue dependence on the input are summarized as 

f Type2 : {ig < ir) D (ig ~ iu) 
\ TypeS : ig < ir 



In Fig.|T3|the above mechanism of the analogue change is shown schematically. 

Note again that our mechanism for input dependence works universally in 
one-way coupled dynamical systems, irrespective of the choice of our specific 
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u 



(b) 



Figure 13: Schematic diagram of the analogue change. Cases with a clear 
analogue change (a), and without it (b). 



model. In our model, the analogue change is represented as the frequency of 
oscillations. This transformation of the input to the frequency can be model 
specific. In general, there may be a variety of ways of transformation of input 
values to dynamics at the down-flow. Still the argument presented here is 
applicable to each case. 



6 Noise effect to input dependence 

In the previous section, the mechanism of input dependence is explained from 
the spatial change of X^{i)- In this section, mechanism of input dependence 
is discussed from a different angle, that is the relation between convective 
instability and noise (fluctuation) intensity. It will be shown that the input 
dependent dynamics at the down- flow is seen only within some range of noise 
intensity. 



Fig.|Tj shows our numerical results on the input dependence of the vari- 
ance of a;*(t) at the down- flow for different amplitudes of noise, where the 
parameters are fixed so that the system belongs to Type3. It is seen that 
the mean square variance has input dependence only in medium intensity (•) 
of the noise. For too large (x) intensity of the noise, the variance does not 
show input dependence, while the dependence gets weaker as noise ampli- 
tude is decreased (+). "Suitable" range of noise intensity is required to have 
analogue change on the input. 
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Figure 14: Mean square variance (< > — < ,x* >2)^/2 is plotted against 
tiie change of input values, for 3 different values of noise amplitudes. + 
((7 = 1.0 X 10-^) ^« ((7 = 3.0 X 10-^) ^x (cr = 8.0 X lO^^). The parameters 
corresponds to Type3, a = 2.7, h = 10.7, c = 4.8, d = 9.4, K = 0.015, e = 4.5. 
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6.1 disappearance of input dependence at a low noise 
regime 

The mechanism of the disappearance of input dependence at a low noise 
regime is straightforwardly explained with the argument of the previous sec- 
tion. According to eg. (]Tl|) , ig gets larger with the decrease of the noise, and 
the relationship between ig and iu or v changes. With the decrease of a, the 
following two changes are possible. 

1. Disappearance of digital change (for Type2): 

With the decrease of a by fixing the input, ig can be larger than i„, and 
the transition from the periodic to stochastic phase, and then to fixed 
point phase occurs. The periodic phase no longer appears irrespective 
of the input values (within their allowed range in our model) for small 
a. Hence digital dependence on inputs disappears. 

Fig.|T3|(a) shows input dependence of the variance at the down-flow 



dynamics with the change of the noise amplitude, while in Fig.|T5|(b) 
input dependence of ig and iu is plotted. It is demonstrated that the 
digital change requires larger input values as the noise amplitude is 
decreased. 



0.03 0.035 0.04 0.045 0.05 0.055 
Input (boundary con 



0.06 0.065 0.07 0.075 



(b) 




Figure 15: Input dependence of the variance at i = 80 (a)and of iu (+) and 
ig (b) are plotted. The noise amplitudes are a = 10~^(x), 10~^(o), 10~^(n), 
10^^(«) for each figure (Note iu does not depend on a). The smaller the noise 
amplitude is, the harder the periodic phase appears within allowed range of 
inputs. Parameters a, b, c, d, K, e are same as in Fig.4. 
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2. Disappearance of analogue change (for Type2 or Type3): 

With the decrease of a, ig increases until it gets larger than v, and the 



change from CaseA to CaseB in § [5 .31 follows. Thus analogue depen- 
dence on the input fades out. 

The lower limit for the noise as to have the input dependence is straight- 
forwardly estimated from the above argument. The noise intensity must 
satisfy 

J ij. (analogue change) , , 

^ \ iu (digital change) ^ ^ 

By introducing the average spatial instability Af = {l/ig)J2T and 
using the expression eq.(|lID for ig, we obtain 



^ ^ ^ _ I exp{—X^ir) (analogue change) 
^ Q^p(^_\Sj^^'^ (digital change) 

6.2 collapse of input dependence for a strong noise in- 
tensity 

The collapse of input- dependence at strong noise is due to a different mech- 
anism. When the noise intensity is too large, the wave, once formed, can 
be destroyed due to the noise at the down-fiowQ The generated oscillation 
is distorted by noise and the transmission to down-flow does not work well. 
In particular, for the Type2 case, due to the convective stability of the fixed 
point at the down- flow, the dynamics can stay near the fixed point, after the 
collapse of the limit cycle. Thus, the input dependence between the periodic 
and fixed point states fades out. There is an upper limit of noise intensity 
aL such that the input dependence disappears for a > a^. 

This upper limit aL cannot be estimated by the convective stability of the 
fixed point. It is a nonlinear effect of noise around the limit cycle oscillation, 
for which no theory is available as yet. 



^The noise here plays the same role as that added to information channel in information 
theory. 
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7 Discussion and conclusions 



In the present paper we have reported boundary (input) dependence in a 
one-way coupled differential equation, and presented a general mechanism 
for it. It is shown that the spatially dependent convective instability and 
noise lead to such boundary dependence. 

Boundary condition dependence has also been studied in dissipative structure||12||. 
However our boundary condition dependence is clearly distinguished from 
earlier studies, and is essential in nature. 

The most important difference is dependence on the system size. The 
boundary condition dependence previously studied is nothing but a finite 
size effect. The dependence disappears as the size gets larger. On the other 
hand, our boundary condition dependence is not due to such a finite size 
effect, but remains even in the infinite size limit. Difference of dynamics, 
created by the boundary (input) difference, is fixed to the down-flow ele- 
ment i and is kept even for the limit i ^ oo. Also it should be noted that 
our boundary dependence is independent of initial conditions, and is stable 
against perturbations. 

Furthermore our mechanism supporting the boundary dependence is novel, 
and should be distinguished from that by the previous studies. In our case, 
it is originated in the spatial amplification of fluctuation by the convective 
instability, in contrast with the use of a stable system in the earlier studies 
1T2| . The amplification rate of fluctuation, characterized by the (local) con- 



vective Lyapunov exponent At,(i), depends on space, and accordingly on the 
boundary (input). After the amplification, the dynamics is stabilized at the 
down-flow. Once stabilized, the dynamics stays as an attractor. Dynamics 
before the achievement of the spatial convergence can depend on the bound- 
ary. On the other hand, the generated difference at the down-flow dynamics 
cannot be changed any more and is fixed, due to the stable dynamics there. 
Hence the boundary dependence is fixed at down-flow. 

We have shown that there are two types of boundary (input) dependence. 
One is the threshold-type (bifurcation-like) change, leading to qualitative 
difference in the down-flow dynamics. In our model, bifurcation from the 
fixed point, to the stable stochastic oscillation, and then to the stable peri- 
odic oscillation is found with the change of the boundary value. The other 
is continuous (analogue) change of down-flow dynamics, depending on the 
boundary value. In our example, the frequency of pulses can continuously 
change with the boundary value. 
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These two types of boundary dependence are quantitatively analyzed by 
introducing three characteristic length scales based on the local convective 
Lyapunov exponent A„(i). The first one is the length scale ig necessary for 
the fluctuation to be amplified to the macroscopic order. The second is the 
length scale iu where the convective instability of the fixed point disappears 
and the dynamics becomes convectively stable. The third is the length scale 
V characterizing the spatial convergence of A^(i). 

For the bifurcation-like change, it is crucial if the amplification of fluctua- 
tion is completed before the fixed point dynamics becomes stable. Otherwise, 
the down-flow dynamics remains at the fixed point motion. Thus the condi- 
tion for the qualitative change against the boundary is determined by ig < i^- 
On the other hand, the analogue change depends on if the amplification is 
completed while the spatial dependence of the convective stability is signifi- 
cant, and thus is judged by the condition ig < i^. Note that these conditions 
also imply that such boundary dependence is seen only within some range of 
noise amplitude a, because ig oc loga from eq.(|TT]). 

Since our mechanism is expressed in universal terms in dynamical sys- 
tems, the input (boundary) dependence is expected to be observed gener- 
ally in a spatially extended dynamical system, as long as the spatially inho- 
mogeneous convective instability exists. It is interesting to search for such 
boundary dependence for open-fluid flow, optical system 0], coupled map 
lattice|Q, ^, and so forth, and to check if the relationship among ig and iu 
and V are satisfied. 

Our motivation for this study is originated in biological signaling prob- 
lems. Then what implications are drawn from our result to such problems? 

In our study the input information is represented by the boundary value, 
while the response is given by a state of the down-flow dynamics. At the 
region i ^ ig required for the wave formation, the input difference is ampli- 
fled by the convective instability. The information on the dynamics at this 
region is translated to the response, through the site dependent Xv{i)- If the 
conditions among ig, iu, and v are satisfled, input information can be trans- 
lated to output information through the chemical dynamics at the signaling 
pathway. 

Now let us come back to the questions raised in §1. The flrst question 
is on the reason for the length and complication of the signaling pathway, 
to which our results have some implication. According to our mechanism, 
we need several sites (larger than ig and iu) to have input dependence. This 
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number of sites depends on the parameter, but cannot be too small, as long 
as the dynamics is not too convectively unstable. Hence the pathway must 
have a sufficient length, to have input dependence. Of course the signal 
pathway does not consist of a single chain of reaction, but several chains are 
mutually influenced. Although extension of our study to a multiple chain 
case remains as a future problem, it is expected that such complication can 
afford an effectively long chain required for the input dependence. 

The second question on the suitable response to inputs is answered by the 
above translation of the input (boundary) to response (down-flow dynamics) 
through spatially dependent convective instability Xv{i)- In particular, the 
digital (bifurcation- like) dependence provides a response with some input 
threshold, while the analogue dependence provides the translation from in- 
put concentration to response frequency. Both the dependence are essential 
to neural and sensory responses. It is also interesting to note that our sig- 
nal transmission mechanism is different from that given by Hodgkin-Huxley 



equation for neural signal transmission . 

The third question on the robust signaling mechanism under thermody- 
namic fluctuation is most clearly answered by our mechanism. Indeed our 
input dependence works in the presence of noise, or rather the noise ampli- 
tude a has to satisfy ai > cr > as- Note that the noise is inevitable in a 
cell system by thermal fluctuation and also due to a relatively small num- 
ber of signaling molecules. Indeed the number of most signal molecules is 
around N ~ 1000. Thus the noise of the order I/a/ZV ~ 0.03 should exist. 
The condition ai > cr > (Js for the input dependence may be suggestive in 
thinking why most cell response works only within some temperature range, 
and works in a system of small number of (signaling) molecules. 

Of course several studies have to be pursued in future, not only as dy- 
namical systems but also for the application of our viewpoint to biological 
signaling phenomena. The following list is under current investigation. 

1. Extension to with complex spatial dynamics of A„(z): In our 
model Xvii) monotonically relaxes to Xv{oo). In general, the relaxation 
can be oscillatory or have chaotic transient, as is found in a one-way 
coupled map lattice (OCML)[|l^. In such case, complex input depen- 
dence may be expected. 

2. Extension to a case with a bi-directional or symmetric coupling: Al- 
though the present uni-directional coupling case provides a most straight- 
forward example to see the convective instability, such instability also 
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exists in a bi-directional coupling case, and is characterized by the co- 
moving Lyapunov exponent A„(z) 0. Indeed some preliminary studies 
show that the boundary condition dependence of our type also exists 
in such case. 

3. Extension to a system with several CS attractors at i — oo. In this case 
choice of an attractor can depend on inputs. Complicated dependence 
on input values may exist in a similar manner with fractal or riddled 
basin structure for a multiple attractor system. 

4. Extension to a non-constant input (boundary). If the boundary os- 
cillates in time, for example, this temporal information of the input 
(boundary) can be translated to response dynamics. Indeed, in a 
OCML|Q, the input oscillation is selectively transmitted to down-flow 
depending on its period. 

5. Extension to a multiple chain system with multiple inputs. In a sig- 
naling process in a cell, several pathways exist, which interact with 
each other. With multiple boundary values (inputs), a new type of 
boundary dependence can be expected, such as a response depending 
on combination of several inputs. 
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A calculation of local co-moving Lyapunov 
exponent of the fixed point {x\^yl) 

Since the fixed points change their value by sites in our model, the con- 
ventional method]^] to calculate co- moving Lyapunov exponent has to be 
extended. To measure the amplification rate of perturbation at one site i, 
we consider a system where {x^,y^) = {xl,yl) for j > i and compute the 
amplification rate per lattice point. Let us consider the evolution equation, 

'^{t) = h^{^{t),^-\t)) (15) 
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with h^{^,^-^) = {f{^,^-^),g{^,^-^)) 
The displacement follows the equation, 

,.^^^OMflls,ii,^?R^Sp-^m (16) 

We solve the equation for ^ = {xl,yl) for j > i and Vt with the initial 
condition 5^(0) = 60 j6. The solution can be written as 

6z^{t) = J^{t)6^{0) (17) 

where 

Mt, . . (18) 

Here, J-'(t) is a 2 x 2 matrix, and is obtained by solving eq.(p!8D under 
the initial condition 5z^{0) = 6oj6. By defining e^i*^*^ and e^2(*) as absolute 
values of the eigen values of J^(t) {A{(t) > A^t)) we obtain, 

XM) = lim - Ai(t) (19) 
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